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. Abstract 

cn 

\^ • We present measurements of the redshift-space three-point correlation function of galaxies in the Sloan 

CO ' Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy 

Op , morphology, color, and luminosity. In order to control systematics due to selection effects, we used r-band, 

volume-limited samples of galaxies, constructed from magnitude-limited SDSS data (14.5 < r < 17.5), and 
further divided the samples into two morphological types (early and late) or two color populations (red 
■ and blue) . The three-point correlation function of SDSS galaxies follows the hierarchical relation well, and 

the reduced three-point amplitudes in redshift-space are almost scale-independent {Qz = 0.5 1.0). In 
Q ■ addition, their dependence on the morphology, color, and luminosity is not statistically significant. Given 

I the robust morphological, color and luminosity dependences of the two-point correlation function, this 

c/3 . implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that a simple 

' deterministic linear relation with the underlying mass could not explain our measurements on these scales. 

Key words: cosmology: large-scale structure of universe — cosmology: observations — methods: 
statistical 
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1. Introduction The determination of the 3PCF of galaxies was pio- 

neered by Peebles and Groth (1975) using the Lick and 
While the first-year WMAP (Wilkinson Microwave Zwicky angular catalogs of galaxies. Groth and Peebles 
Anisotropy Probe) data imply that the primordial density (1977) found that the 3PCF C(^i2, ''23, ''31) obeys the fol- 
fluctuations at z « 1000 are virtually Gaussian (Komatsu lowing hierarchical relation: 



C(''l2,f23,f3l) 

= Qr K(r-12)e(r-23) + e(r23)e('^3l) +^(r3l)e('^12)], (1) 



et al. 2003), the present structure traced by galaxies 
shows significant non-Gaussianity. Such non-Gaussian 
features naturally arise during the course of nonlinear 

gravitational evolution of dark-matter density fields, the with being a constant, and ^(r) is the two-point cor- 

formation of luminous galaxies and their subsequent evo- relation function (2PCF). The value of Qr in real space 

lution. Therefore, the degree and nature of non-Gaussian de-projected from these angular catalogues is 1.29 ±0.21 

signatures in the galaxy distribution will provide impor- for r ^ 3/i~^Mpc. Although subsequent analyses of red- 

tant empirical constraints on the physics of galaxy forma- shift catalogs confirmed the hierarchical relation, at least 

tion. The three-point correlation function (3PCF) is the approximately, the value of Qz (in redshift space) appears 

lowest-order unambiguous statistic to characterize such to be smaller, 0.5 ^ 1 (Bean et al. 1983; Hale-Sutton et 

non-Gaussianities. al. 1989; Jing, Borner 1998). Recently, Jing and Borner 

Z Z I Ta J- t- c 'D\, ■ jaj-u- (2003) have studied the luminosity dependence of Q, and 

Present address: Department oi Physics and Astrophysics, ^ ' j f ' 

Nagoya University, Nagoya 464-8602 found a Small, but significant, trend that brighter galaxies 

tend to have a lower amplitude than fainter ones. 
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N-body simulations and perturbative analysis generally 
predict that the 3PCF should depart from the hierarchical 
relation, especially on non-linear scales (Fry 1984; Suto, 
Matsubara 1994; Jing, Borncr 1997; Barriga. Gaztafiaga 
2002). However, Matsubara & Suto (1994) demonstrated 
that redshift distortions substantially reduce the scale- 
dependence of Qr, resulting in Qz {Q in redshift space) 
being nearly constant. Still, the amplitude of the Qz of 
galaxies is roughly (50% 100%) smaller than predicted 
by N-body simulations. The discrepancy is likely to be 
the effect of biasing. Recently, Takada and Jain (2003) 
proposed a phenomenological model to predict the 3PCF 
as a function of the galaxy properties, which is based on 
the halo model. They argue that the color dependence 
of Qr should be strong, i.e., the galaxy biasing severely 
affects the 3PCFs. 

There are several statistics that are closely related to 
the 3PCF, e.g., the reduced moments of counts in cells, 
such as skewness and kurtosis (see Szapudi ct al. 2002, 
2004 in preparation). Another similar statistics is the bis- 
pectrum, which is the Fourier counterpart of the 3PCF. 
Recently, Verde et al. (2002) computed the bispectrum of 
the 2dF Galaxy Redshift Survey, and concluded that the 
non-linear bias is consistent with zero. The other com- 
plementary approach to quantify the non-Gaussianity in 
the SDSS galaxies by means of topological analysis (i.e., 
in terms of the Minkowski functionals) was already con- 
ducted (Hoyle et al. 2002; Hikage et al. 2002, 2003). The 
above results are largely consistent with that predicted 
from purely gravitational nonlinearity, while a weak mor- 
phological difference is marginally detected (e.g., figure 13 
of Hikage et al. 2003). Here, we report on the first results 
of characterizing the non-Gaussianity in the SDSS galaxy 
samples using 3PCFs in redshift space. In particular, we 
consider the dependence of 3PCFs on the morphology, 
color, and luminosity of galaxies by constructing volume- 
limited samples so as to avoid possible systematics due to 
the selection function. A separate companion paper by 
Nichol et al. (in preparation) focuses on a detailed com- 
parison of 3PCFs between the 2dF (Jing, Borner 2004) 
and the SDSS galaxies. 

The paper is organized as follows: section 2 describes 
our volume-limited sample of SDSS galaxies and the mor- 
phology and color classification methods. The measure- 
ments of 3PCFs are detailed in section 3 with particular 
attention given to their dependence on the morphology, 
color and luminosity of galaxies. Finally, we summarize 
the results in section 4. We compare several different es- 
timators of 2PCFs and 3PCFs in Appendix. Throughout 
the data analysis, we adopt the following set of cosmolog- 
ical parameters: the matter density parameter, Jim = 0.3; 
the dimensionless cosmological constant, fiA = 0.7; and 
the Hubble constant in units of lOOkm/s/Mpc, h ~ 0.7. 

2. Volume-Limited Samples of SDSS Galaxies 

It is expected that galaxy clustering depends on the 
intrinsic properties of the galaxy samples under consid- 
eration, including their morphological types, colors, and 



luminosities. Therefore, a straightforward analysis of 
magnitude-limited samples must be interpreted with cau- 
tion because several different effects may be simultane- 
ously involved. The magnitude-limited sample of SDSS 
galaxies enables us to construct volume-limited samples 
of different luminosities, morphologies and colors. This 
essentially removes the above difficulties while keeping a 
statistically significant number of galaxies in each sample. 

Our present analysis is based on a subset of the SDSS 
galaxy redshift data, 'Large-scale Structure Sample 12' 
(Blanton et al. in preparation), which is larger by a fac- 
tor of 1.8 than the public Data Release One (Abazajian 
et al. 2003). This sample includes galaxies with r-band 
magnitudes of between 14.5 and 17.77 after correction for 
Galactic reddening using the maps of Schlegel, Finkbeiner, 
and Davis (1998). For more details, the following can 
be consulted: York et al. (2000) for an overview of the 
SDSS; Gunn et al. (1998) for a description of the pho- 
tometric camera; Stoughton et al. (2002) for photometric 
analysis; Fukugita et al. (1996), Hogg et al. (2001), and 
Smith et al. (2002) for the photometric system; Pier et 
al. (2003) for the astrometric calibration; Eisenstein et 
al. (2001); Strauss et al. (2002) for selection of the galaxy 
spectroscopic samples; and Blanton et al. (2003) for spec- 
troscopic tiling. 
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Fig. 1. Redshift - absolute magnitude relation for vol- 
ume-limited samples of galaxies with 14.5 < r < 17.0. 

We construct two different kinds of volume- limited sam- 
ples; one is based on the morphological classification, and 
the other is on the color of galaxies (figure 1 and table 1). 
Figure 2 illustrates the distribution of those galaxies in a 
volume-limited sample. 

Following Shimasaku et al.(2001), we classify the mor- 
phology of galaxies according to the (inverse) concen- 
tration index (q), which denotes the ratio of the half- 
light Petrosian radius to the 90%-light Petrosian radius. 
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Fig. 2. Distribution of galaxies in tlic volume-limited sam- Fig. 3. Completeness of morphological classification ac- 

ples (14.5 < r < 17.0; see table 1 and figure 1). cording to the inverse concentration parameter. 



Table 1. Volume-limited samples of 
14.5 < r < 17.0 for the early/late 
14.5 < r < 17.5 for color classification. 



SDSS galaxies; 
classification and 
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red 


blue 


-22 to - 
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0.065-0.14 
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8329 


-21 to - 


-20 


0.042-0.093 


8930 
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-20 to - 
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0.027-0.061 
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Specifically, we define galaxies with a < ci^c as early-types, 
and those with Ci > Ci^c as late-types. The critical value 
Ci^c = 0.35 is adopted for galaxies with r < 16.0. Figure 3 
shows the completeness and contamination of this classifi- 
cation scheme. Around Ci 0.35 the completeness reaches 
80% and the contamination goes below 20%. The clas- 
sification based on Ci is affected by the fact that the 
image quality of SDSS galaxies with r > 16.0 starts to 
be degraded due to the seeing. In order to empirically 
compensate for the effect, we slightly change the critical 
value: a^c 0.359 for 16.0 <r < 16.5 and a^c = 0.372 
for 16.5 < r < 17.0. Those values are chosen so that the 
number ratio of early- to late-types remains the same in 
the three different magnitude ranges: 14.0 < r < 16.0, 
16.0 < r < 16.5, and 16.5 < r < 17.0. We do not classify 
the morphology of galaxies with r > 17.0. 

On the other hand, classification of galaxies according 
to their colors does not require such a good image quality 
as the Ci classification. Thus, we can make use of galax- 
ies down to r — 17.5, leading to an improved statistical 
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Fig. 4. Histgrams of the color distribution for three vol- 
ume-limited galaxy samples (table 1). The bin of his- 
togram is A(g — r) = 0.025. The dotted vertical lines 
in each panel indicate the color thresholds that di- 
vide galaxies in three different color subsamples. 

significance. Figure 4 shows the distribution of g — r for 
volume-limited galaxies corresponding to the lower part of 
table 1 . Clearly, the colors are strongly correlated with the 
luminosity, and several different versions of color-selected 
samples are possible. Since our goal is to sec if galaxy 
clustering is dependent on the colors in a complementary 
manner to the previous morphological classification, we 
divided each volume-limited sample into three subsam- 
ples that roughly have equal numbers of galaxies. More 



4 



I. Kayo et al. 



precisely, we defined red (blue) galaxies as g — r > 0.78, 
> 0.75, and > 0.72 {g - r < 0.69, < 0.63, and < 0.55) for 
-22 < Mr ~h\ogh < -21, -21 < Mr - 5\ogh < -20, and 
—20 < Mr — 5\ogh < —19, respectively. Here, g and r are 
k-corrected magnitudes. These thresholds are indicated 
in the vertical lines in figure 4. Figure 5 illustrates the 
relations between our two classification schemes. 
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Fig. 5. Relation between the colors and the inverse con- 
centration parameters of galaxies. The horizontal dotted 
lines correspond to the vertical dotted lines in figure 4. 
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Fig. 6. Parameters to define the shape of triangles. 



3. Measurements and Results 

3.1. Counting Triplets 

In the present analysis, we directly computed the num- 
ber of triplets to estimate the 3PCFs C(si2, S23, S31) in 
redshift space. In this way, we properly take account of 
the complicated survey volume boundary. To be more spe- 
cific, we distributed random particles over the survey 
volume of a given sample of Njj galaxies, following the 
selection function of the latter (we used A^r = 25iVD.aii, 
where A^D.aii is the number of all galaxies in each lumi- 
nosity bin). Then, the 3PCFs were computed from the 
following estimator (Szapudi, Szalay 1998): 

C(si2,S23,S3l) 

RRR{si2,S23,S3i) RRR{si2,S23,S3l) 
^ ^ DRR{si2,S23,S3i) ^ 

RRR{si2,S23,S3l) 

where DDD, DDR, DRR, and RRR are the number 
of corresponding triplets consisting of galaxies and ran- 
dom particles, and are normalized by Nd{Nd — 1){Nd — 
2)/6, Nd{Nd~1)Nr/2, NDNR{Nn-l)/2, andNiiiNn- 
^){Nr — 2)/6, respectively. The above estimator is close 
to optimal, because it is constructed to provide accurate 
edge effect corrections. We decided to use the Szapudi and 
Szalay estimator, since we found that it converges the true 
3PCF slightly more rapidly than the estimator of Jing and 
Borner (1998) with the given number of random particles 
(see Appendix for details). We also estimated the redshift 
space 2PCF with the analogous estimator (Landy, Szalay 
1993), and defined the normalized amplitude of 3PCFs as 

Qzisi2,S23,S3l) 

^ Cisi2,S23,S3l) 

~ e(si2)?(s23) +e(s23)e(s3l) +?(s3l)e(si2) ' ^ ' 

Two different parameterizations of triangular shape are 
conventionally used in this field. One is 

(4) 

where S12 < S23 < S31 is assumed (thus m > 1 and < w < 1). 
The other is 

S12 = s, S23 = qs, S31 = sy/l + q'^ - 2qcos0, (5) 

where S12 < S23 is assumed (thus, g > 1, but S31 may be 
smaller than either of the other two; see figure 6). Note 
that according to the second parameterization, one trian- 
gle is counted in three different bins in {s,q,6). Although 
the first parameterization is more traditional in the real 
data analysis, we adopt the second, which is widely used 
in theoretical predictions. More specifically, we choose 
8 equally-spaced logarithmic bins in 0.4/i-^Mpc < s < 
10.0/i^^Mpc and 5 equally-spaced linear bins both in 
1 < g < 5 and in < < tt. Since the number of triplets 
of s < 1.0/i~^Mpc is very small (typically DDD ^ 50), we 
plot the 3PCFs only for s > 1.0/i~"'^Mpc (where we have 
more than 100 DDD counts) in what follows. Also, the 
errors quoted below are evaluated from the 16 jack-knife 
re-sampling (sec e.g., Lupton 1993). 




S12 = s, S23 = us, S31 = {u + v)s, 
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3.2. Equilateral Triangles 

We ignore any possible dependence on the triangular 
shape for the moment, and first consider 3PCFs for equi- 
lateral triangles in detail. For this purpose, we first di- 
rectly evaluate the lengths of three sides of triplets. If 
those sides fall in the same bin among 8 equally-spaced 
logarithmic bins between 0.4ft,~^Mpc and 10.0ft.~^Mpc, we 
define the triplets as being equilateral triangles. Figures 
7, 8, and 9 show over l/i^^Mpc <s< lO/i^^Mpc. The 
overall conclusion is that Qz is almost scale-independent 
and ranges between 0.5 and 1.0, and no systematic depen- 
dence is noticeable on the luminosity, morphology, and 
color. 

Previous simulations and theoretical models (Suto 
1993; Matsubara, Suto 1994, Matsubara 1994, Takada, 
Jain 2003) indicate that Q decreases with the scale in 
both real and redshift spaces. Observationally, Jing and 
Borner (2003) suggested that Qz of the 2dF galaxies also 
decreases as scale. Neither trend is clear in the results. 
This might be partly due to the different velocity disper- 
sion of the galaxies in the two samples. Jing, Mo, and 
Borner (1998) reported that the early spirals (Sa, Sab, 
Sb) have a significantly smaller Q than the other types, 
which is qualitatively consistent with the present results. 
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Fig. 8. Same as figure 7, but different symbols corre- 
sponds to different morphologies; all galaxies in open circles, 
early-type galaxies in solid triangles, and late-type galax- 
ies in crosses. Top: -22 < Mr - 51og/i < -21, Middle: 
-21 < Mr - 51ogh < -20, Bottom: -20 < - 51ogh < -19. 
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Fig. 7. Luminosity dependence of the normalized ampli- 
tude of three-point correlation functions classified accord- 
ing to the morphologies of galaxies. Only equilateral tri- 
angles are considered. Different symbols correspond to dif- 
ferent luminosities: —22 < M,. — 51ogh < —21 by solid tri- 
angles, —21 < M-i 5 log h < — 20 by open circles, and 

—20 < Mr — 51og/i < —19 by crosses. Top: all galaxies. 
Middle: early-type galaxies, Bottom: late-type galaxies. 
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Fig. 9. Same as figure 8, but for classification according 
to the colors of galaxies. All galaxies by open circles, red 
galaxies by solid triangles, and blue galaxies by crosses. 



3.3. Relation to Biasing of the Two-Point Correlation 
Function 

The statement that the values of Qz are insensitive to 
the intrinsic properties of galaxies may sound somewhat 
misleading. The reduced three-point correlation ampli- 
tude is normalized by the square of the 2PCFs for the 
corresponding galaxies. This implies that the (unreduced) 
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3PCF depends on the galaxy properties, since 2PCFs 
are known to show a clear dependence on the luminos- 
ity, color, and morphology of galaxies (e.g., Lahav and 
Suto 2003 for a recent review) . In order to further demon- 
strate the expected dependence in the current samples , 
we compute the biasing parameters estimated from the 
2PCFs, 



Cz,acdm(s) ■ 



(6) 



where the index i runs over each sample of galaxies with 
different morphologies, colors, and luminosities. The pre- 
dictions of the mass 2PCFs in redshift space, C2,acdm(s), 
in the A Cold Dark Matter model are computed by prop- 
erly taking account of the light-cone effect over the corre- 
sponding redshift range (Hamana et al. 2001). In doing 
so, we adopt the fluctuation amplitude as = 0.9 and the 
one-dimensional peculiar pairwise velocity dispersion of 
800km/s in addition to our fiducial set of cosmological 
parameters. 

As an illustrative example, consider a simple bias model 
in which the galaxy density field, 6g,i, for the i-th popu- 
lation of galaxies is given by 

'^g,' = bg.i{l)Sniass + feg,i(2)'5mass- (7) 

If both 6g^i(i) and &g,i(i) are constant and the mass density 
field (5mass ^ 1, equation (3) implies that 



i(l) 



62 



(8) 



Thus, the linear bias model (6g,i(2) = 0) simply implies 
that Qg.i is inversely proportional to b^ ^iy While this 
simple model may not be accurately applicable on the 
scales of our results, it is instructive to plot l/6g,i de- 
duced from the 2PCF of galaxies with different lumi- 
nosities and morphologies. The results are shown in fig- 
ures 10, 11, and 12. A linear bias model would pre- 
dict that the ratio of of early-types and late-types 
is Qz,cariy/Qz,iatc = ^latc/^oariy ~ 0.8, and similarly that 
(3z,rcd/Qz,biuo = buuc/brcd ~ 0.5. Neither figure 8 nor 9, 
however, shows such systematic trends within our 20 
percent measurement accuracy. In a sense, the biasing 
in the 3PCFs seems to compensate for the difference of 
Qg purely due to that in the 2PCFs. Such a behavior is 
unlikely to be explained by any simple model inspired by 
the perturbative expansion, like equation (7). Rather, it 
indeed points to a kind of regularity or universality of the 
clustering hierarchy behind galaxy formation and evolu- 
tion processes. At least we can conclude that the galaxy 
biasing is more complex than the simple deterministic and 
linear model. More precise measurements of 3PCFs, and 
even higher order statistics with future SDSS datasets 
would indeed be valuable to gain more specific insights 
into the empirical biasing model. 

3.4- Shape Dependence 

Next consider the dependence of the 3PCFs on the 
triangular shape. For this purpose, we consider the 
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Fig. 10. Inverse of the biasing parameters of tiic two-point 
correlation functions plotted in the same way as in figure 7. 
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Fig. 11. Inverse of biasing parameters of the two-point cor- 
relation functions plotted in the same way as in figure 8. 



volume-limited samples classified according to the colors 
(r < 17.5) because they have a greater number of galaxies. 
Figure 13 shows Qz for all, red, and blue galaxies with 
—21 < Mr — 51og/i < —20 separately. Those plots indi- 
cate a weak ^-dependence expected from the previous per- 
turbation theory and N-body simulations (e.g., Barriga, 
Gaztanaga 2002; Takada, Jain 2003). The 6'-dcpendence 
is weaker on nonlinear scales (s < 1/i^^Mpc), and becomes 
noticeable on larger scales. While the amplitudes of Qz 
are much smaller than, and the color dependence seems 
very different from the corresponding prediction (figure 
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Fig. 12. Inverse of biasing parameters of the two-point cor- 
relation functions plotted in the same way as in figure 9. 

17 of Takada, Jain 2003), they adopt a simple halo ap- 
proach, and our color selection criteria are not necessar- 
ily the same as theirs; further careful study is needed to 
understand the origin of the difference (after the present 
paper is accepted, Wang et al. 2004 submitted a theo- 
retical paper applying the halo model to account for the 
3PCFs of 2dF galaxies, in particular their dependence on 
types and luminosity dependences, and shows the possi- 
bility that the halo model can match the observation). 
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Fig. 13. Shape dependence of Qz of vol- 
ume-limited samples of galaxies with r < 17.5 and 
—21 < Mr — 51og/i < —20 (classified according to their color). 



4. Summary and Discussion 

We have presented the first detailed study of the three- 
point correlation function of SDSS galaxies in redshift 
space. We examined the dependence of their reduced am- 
plitudes, Qz, [eq.(3)] on the scale and shape of the tri- 
angles for various volumed-limited samples with different 
morphologies, colors, and luminosities. As for equilateral 
triangles, wc basically confirm the hierarchical clustering 
relation for scales of l/i~^Mpc < s < 10/i~^Mpc. On the 
other hand, Qz shows significant shape dependence, as ex- 
pected from perturbation theory and N-body simulations, 
particularly at large scale, although the amplitude of Qz 
is much smaller than the expectation. 

The most important finding is that although the two- 
point correlation functions of those galaxies show a clear 
dependence on their luminosity, morphology, and color, 
we did not find any robust dependence of Qz on those 
properties; the Qz results when taken alone imply that 
galaxies are faithful tracers of the underlying mass distri- 
bution, and that their intrinsic properties are independent 
of the purely nonlinear gravitational evolution. At the 
same time, however, their two-point correlation functions 
strongly indicate that they are biased tracers. If we at- 
tempt to reconcile these two behaviors by using a simple 
linear bias model, we cannot adjust it to simultaneously 
explain both the two-point correlation functions and the 
insensitivity of Qz on the galaxy properties. This implies 
that the galaxy biasing is non-linear and fairly complex, 
while nonetheless displaying a remarkable constancy in 
reduced amplitudes. 

The above conclusion may appear to be inconsistent 
with the argument of Verde ct al. (2002), who concluded 
that the non-linearity in 2dF galaxy biasing is negligi- 
ble from the bispcctrum analysis. We note, however, that 
they did not attempt to divide those galaxies into different 
subsamples. If we look at the results of all galaxies alone, 
top-panels of figures 7 and 10 may indicate the same con- 
clusion as reached by Verde et al. (2002). In other words, 
our empirical finding indicates that the non-linearity of 
the biasing, &2, of galaxies vanishes when the correspond- 
ing linear biasing parameter, &i, is unity. The fact that 
the two parameters, bi and 62, are highly correlated with 
each other is qualitatively consistent with some biasing 
model predictions including halo models or peak models 
(Mo et al. 1997). It is remarkable that the current galaxy 
survey approaches the level of discussing the validity of 
those models quantitatively. Future SDSS datascts with 
more area and a larger number of galaxies will be able 
to place stronger limits on this observation. This in turn 
will be useful for a detailed understanding of the nature 
of galaxy biasing. 
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Appendix. Comparison of the Estimators of 
Two- and Three-Point Correlation Functions 

For an accurate determination of the normalized am- 
plitudes of Q, we need reliable estimators for the two- 
and three-point correlation functions. For this purpose, 
we use a set of data particles (A^^i = 64^; generated 
from a N-body simulation with 256'^ dark matter parti- 
cles; Jing and Suto 1998) distributed in a cube of Vbox — 
(100/i~^Mpc)^. We then compare the following three es- 
timators for 2PCFs: 



DD 



direct 



1, (direct); 



(Al) 



1^2 /T4c 

^«"i§i?F"^' (Hamilton 1993); (A2) 
DD-2DR + RR 



RR 



, (Landy, Szalay 1993); (A3) 



where V12 denotes the volume of the spherical shell with 
the thickness of the corresponding separation bin, and 
DD, DR, and RR are the number of corresponding pairs 
of data and random particles, and are normalized by 
Nd{Nd - l)/2, NdNr, and Nr{Nr - l)/2, respectively. 

Similarly, we use the following estimators for the 
3PCFs: 



C,IB = 







_ DDD 

^123/ ^box 

[RRRfDDD 3RRRDDR 
[DRR]^ [DRRp ^ ' 

DDD - 3DDR + 3DRR - RRR 



ss 



[Vol. , 
(A4) 

(A5) 
(A6) 



RRR 

The first uses the analytical expression for the effective 
volume squared for the equilateral triplets. If one consid- 
ers the size of the the equilateral triplets between rmin and 
fmax, it is simply given by 



2/„2 
max 



(A7) 



Because similar expressions in an arbitrary survey volume 
shape are not available analytically, we focus on the equi- 
lateral triplets in periodic cube of simulations. The latter 
two correspond to estimators by Jing and Borner (1998) 
and Szapudi and Szalay (1998), respectively. The com- 
bination of terms in Jing and Borner (1998) comes from 
the fact that the first term corresponds to 1-1-^(^12) + 
^('"23) +^(r3i) -|-C(ri2,r23,r3i), and the second term to 
3 + ^(^'12) + ^(^23) + C(r3i), yielding the desired ( mathe- 
matically. 

The results of the convergence test are shown in Figure 
14 for 2PCFs and Figure 15 for 3PCFs. These figures show 
the ratios of the two different estimators with various 
against the "direct" estimates [equations (Al) and (A4)] 
which do not use the random particles. With this number 
of random particles, equations (A3) and (A6) reproduce 
the "direct" results within ~ 5% for 2PCFs and 3PCFs, 
respectively. 

Figure 14 indicates that ^ls converges to the true value 
(^direct) morc rapidly even with a smaller Nn than 
(see also Kerscher et al. 2000). Figure 15 shows that the 
estimators of 3PCFs are generally less stable compared 
with those of 2PCFs. Still, Css seems better behaved than 
Cjb. Thus, we adopted ^ls and ^ss in the present analysis. 
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